Charge pumping in magnetic tunnel junctions: Scattering theory 
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We study theoretically the charge transport pumped by magnetization dynamics through epitaxial 
FIF and FNIF magnetic tunnel junctions (F: Ferromagnet, I: Insulator, N: Normal metal). We 
predict a small but measurable DC pumping voltage under ferromagnetic resonance conditions for 
coUinear magnetization configurations, which may change sign as function of barrier parameters. A 
much larger AC pumping voltage is expected when the magnetizations are at right angles. Quantum 
size efi'ects are predicted for an FNIF structure as a function of the normal layer thickness. 



A magnetic tunnel junction (MTJ) consists of a thin 
insulating tunnel barrier (I) that separates two ferromag- 
netic conducting layers (F) with variable magnetization 
direction.^ With a thin normal metal layer inserted next 
to the barrier, the MTJ is the only magnetoelectronic 
structure in which quantum size effects on electron trans- 
port have been detected experimentally^. More impor- 
tantly, MTJs based on transition metal alloys and epi- 
taxial MgO barriers^ are the core elements of the mag- 
netic random-access memory (MRAM) devices^ that are 
operated by the current-induced spin-transfer torque.^:^. 

It is known that a moving magnetization of a ferromag- 
net pumps a spin current into an attached conductor!^ 
Spin pumping can be observed indirectly as increased 
broadening of ferromagnetic resonance (FMR) spectrai^ 
The spin accumulation created by spin pumping can be 
converted into a voltage signal by an analyzing ferro- 
magnetic contact. This process can be divided into two 
steps: (1) the dynamical magnetization pumps out a spin 
current with zero net charge current, (2) the static mag- 
netization (of the analyzing layer) filters the pumped spin 
current and gives a charge current. In the presence of 
spin-flip scattering, the spin-pumping magnet can gener- 
ate a voltage even in an FN bilayerJii^ Spin-pumping 
by a time-dependent bulk magnetization texture such as 
a moving domain wall is also transformed into an elec- 
tromotive force.^.'^ Other experiments on spin-pumping 
induced voltages have also been reported i^'*'^^ 

Here we present a model study of spin-pumping in- 
duced voltages (charge pumping) in MTJs. Since the 
ferromagnets are separated by tunnel barrier, we can- 
not use the semiclassical approximations appropriate for 
metallic structures jiSiiiii^iii Instead, we present a full 
quantum mechanical treatment of the currents in the 
tunnel barrier by scattering theory. The high quality of 
MgO tunnel junctions and the prominence of quantum 
oscillations observed in FNIF structures (even for alu- 
mina barriers) provide the motivation to concentrate on 
ballistic structures in which the transverse Bloch vector 
is conserved during transport. For a typical MTJ un- 
der FMR with cone angle 9 ~ 5° at frequency / = 20 
GHz, wc find a DG pumping voltage of |Vcp| — 20 nV 
for coUinear magnetization configurations or AG voltage 
with amplitude Vcp — 0.25 fiY for perpendicular config- 
urations. The magnetization dynamics-induced voltages 



could give simple and direct access to transport parame- 
ters of high-quality MTJs, such as barrier height, magne- 
tization anisotropics and damping parameters in a non- 
destructive way. The polarity of the pumping voltage 
can be changed by engineering the device parameters, 
etc. An oscillating signal as a function of the thickness 
of the N spacer leads to Fermi surface calipers that are in 
tunnel junctions not accessible via the exchange coupling. 

We consider a structure shown in Fig.[lja), where two 
semi-infinite F leads (F(L) and F(R)) are connected by 
an insulating layer (I) of width d and a non-magnetic 
metal layer (N) of width a. The magnetization direc- 
tion ofF(L)/F(R), mi/m2 (|mi| = |m2| = 1), is treated 
as fixed/free. We disregard any spin accumulation in F, 
thus treat them as ideal reservoirs in thermal equilibrium. 
This is allowed when the spin pumping current is much 
smaller than the spin- flip rate in the ferromagnet, which 
is usually a good approximations. The structure reduces 
to an FIF MTJ when a = 0. Let A, B, . . . , F he the spin- 
dependent amplitudes {A^ — (^|,yl|)) at specific points 
(see Fig. [1]) of /lux-normafeed spinor wave-functions. The 
scattering states can be expressed in terms of the incom- 
ing waves A and F, such as: 

E = seaA + sefF, (1) 
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FIG. 1; (Color online) (a): FNIF heterostructure, in which 
Si, 2 indicate two different scattering regions, (b): Potential 
profiles for majority and minority spins in F are shown by 
solid and dashed lines. The exchange splitting is A and the 
tunnel barrier has height 17b relative to the Fermi energy Ey- 



where sea and sef are 2x2 matrices in spin space and 
can be calculated by concatenating the scattering matri- 
ces of region 51,2 and of the bulk layer I. To first order 
of the transmission (ib) through the bulk I, 



SEA = i2 [(1 - r;f2)-4b(l - r'lrb)-^] h, 



(2) 



where ii, 2/^^1,2 are the 2x2 transmission/reflection ma- 
trices for Si^2 (see Fig.[T]), the hatless tb/^b are the spin- 
independent transmission/reflection coefficient for the in- 
sulating bulk I. The primed and unprimed version spec- 
ify the scattering of electrons emitted coming from the 
left and right, respectively. The reflection coefficient rb 
is due to the impurity scattering inside the bulk I, and 
its magnitude mainly depends on the impurity density 
in I, especially near the interfaces. All scattering coeffi- 
cients are matrices in the space of transport channels at 
the Fermi energy that are labeled by the transverse wave 
vectors in the leads: q, q' (the band index is suppressed). 

The response to a small applied bias voltage can be 
written as Jc ^ GcV with conductance Gc'. 
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where Tro- [• • •] denotes the spin trace and the summation 
is over all transverse modes in the leads at the Fermi level. 

When the structure is unbiased but the magnetic con- 
figuration is time-dependent, a spin current is pumped 
through the structure?- When the dynamics is slow, 
rhi ^ Ep I h, it can be treated by the theory of adiabatic 
quantum pumpingji^ We consider a situation in which 
the magnetization (1x12) of one layer precesses with ve- 
locity (j) around the z-axis with constant cone angle 9, 
whereas the other magnetization (mi) is constant (see 
Fig. [T|) . We focus on the charge current that accompa- 
nies the spin pumping: 
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When a DC current is blocked (open circuit), a voltage 
bias Vcp builds up 
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The discussion above is valid for general scattering 
matrices that e. g. include bulk and interface disor- 
ders. In order to derive analytical results, we shall make 
some approximations. First of all, we assume that spin 
is conserved during the scattering, U for Si {i = 1,2, 
similar for fi) is collinear with 111^:^ Expanded in Pauli 
matrices a = (0-3;, (jy, ct^), ii ^ + t^o- ■ nij, with 
tf = {tj ±t\)/2. tf (a =1,1) is the transmission ampli- 
tude for spin-up/down electrons with spin quantization 



axes xiii in the scattering region Si. In the absence of 
impurities (rb = rf^ = 0), Eq. ^ becomes 

SEA = {tt^htt + ^^'otimi ■ m2) (6) 
+ (T • (<Jib<]^mi + <2^tbi]'"m2 — ztj^ib^r™! ^ "^2) ■ 

Since all hatless quantities in this equation are still ma- 
trices in k-space, such as = tj(q, q')j the order of 
i2,ib,ii as in Eq. ^ should be maintained. The sef 
term in Eq. ^ may be disregarded, because only the 
part of Sep that depends on both mi and m2 contributes 
to jcp, and that part is in higher order of tb. 

Another approximation is the free electron approxima- 
tion tailored for transition metal based ferromagnetsiii^ 
We assume spherical Fermi surfaces for spin-up and spin- 
down electrons (in both F(L) and F(R)) with Fermi wave- 
vectors kl = ^2mEYlh^ and k\, = ^2m{E^ - 
with an effective electron mass m in F. Electrons in N are 
assumed to be ideally matched with the majority elec- 
trons in F (fcp = ^f,wn = "rn). Let C/b and mb = /3m 
be the barrier height of and effective mass in the tun- 
nel barrier. The adopted potential profile is shown in 
Fig. [T]^b). We assume the transverse wave- vector q to 
be conserved (q = q') by disregarding any impurity or 
interface roughness scattering, which means the scatter- 
ing matrices (i? 21^1^2; ^b) are diagonal in k-space. With 
these approximations, the double summation in Eqs. ([Sj 
3]) is replaced by a single integration over transverse wave- 
vectors. The scattering amplitudes and rf can be cal- 
culated by matching the flux-normalized wave-functions 
at the interfaces. The transmission coefficient in the 
barrier bulk is the exponential decay: tb — e~'^'^ with 
K = yj2/m\Jj\jW~+l^ . Then we obtain our main result 
from Eq. ([6]): 



3c = ^e"'"' {T+T+ + Tf T2-mi • m2) , 
Jcp = ^e-^-'' T-mr 

[|t^P(m2 xm2)+Im(4*<^)m2] , 
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where 7;+ = |tj|2 + |ij|2 is the total tr ansmission probabil- 



ity for scattering region Si , and = pi = 1 4 P ~ I ^ 
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with polarization pi = T~ /T+ . In Eq. jTbl, The term 
in the square brakets is the transmitted spin pumping 
current, and Tfmi represents the filtering by the static 
layer that converts the spin into a charge current. 

For an Fe/MgO/Fe MTJ: fcj, = 1.09 k'^ and 4 = 0.42 
for Fer^ and J7b — 1 eV and (3 = m^/m = 0.4 
for MgO,^ This implies Ey ~ 4.5 eV, A ~ 3.8 eV 
« 0.85£;f, and Ub « 0.25SF (<b < 1 when Uy, > O.lEp 
and c? > 0.5 nm). For an FIF structure (a = 0), both Si 
and 5*2 contain only a single F(L)/I (for Si) or I/F(R) 
(for S2) interface. From the potential profile in Fig.[T]^b), 
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barrier height: UJE^ barrier height: UJE^ insertion layer width: a (nm) 



FIG. 2: (Color online) Barrier height (Uh) and width (d) de- FIG. 3: (Color onUne) Vcp vs. N layer thickness a for FNIF 
pendence of the pumping voltage. Left: DC. Right; Maxi- ([/^ ~ 0.25£'f). 
mum amplitude of AC voltage. 



for fcj^ — /cp^ — (7^ > and zero otherwise. 

If m.2 precesses about an axis that is parallel to mi 
(x = 0° or 180°, see Fig. [IJ mi • 012 = and the second 
term in Eq. (j7bp vanishes. The dot product |mi • {m.2 x 
m2)| = 27r/sin^0 is time-independent, thus generates a 
DC signal. Let us consider an FIF MTJ with barrier 
width d = 0.8 nm, with m2 precessing around the z-axis 
at frequency f ~ 20 GHz with cone angle 9 = 5°. We 
find a DC charge pumping voltage over the F leads V^p ~ 
15 nV when mi is parallel to the precession axis (x = 
0°) and Vcp ~ -19 nV when anti-parallel (x = 180°).2i 
\Vcp\ is higher for the anti-parallel configuration simply 
because its resistance is higher. When the precession 
cone angle 6* = 10°^ the DC voltage |Vcp| ~ 60 nV, 
similar to a previously measured pumping voltage in a 
metallic junctioniiS 

Fig.dja) shows the DC Vcp as a function of the barrier 
height Uh for an FIF structure a.t 6 = 5° and / = 20 
GHz. Vcp = Gc^Jcp increases as a function of barrier 
height mainly because Uc^jcp increases as a function of 
Uw- From Eq. Q, we have g-^jcp « [T^ lT+){\t2? IT+) 
(assume T{T^ < T+T+). The first ratio Tf /T+ = 
Pi oc [[h/PY — + ^T^J' increases as a 

function of Uh through K.{Uh), whereas second ratio 
It^P/Tj^ oc {yjk^ — y/klY/lk-^ + is independent of 
The pumping voltage therefore increases with Uh 
(and 1//3). We also see that Vcp decreases when d in- 
creases, which can be understood by the following: The 
effect of the tunnel barrier is to focus the transmission 
electrons on small q's due to the exponential decay fac- 
tor exp[— 2K(q)(i]. Smaller q implies larger kinetic energy 
normal to the barrier and therefore reduced sensitivity to 
the spin-dependent potentials. Hence, Vcp decreases with 
barrier width. The lowest curve in Fig. [DJa) is approx- 
imately Vcp = 5(r^('Z)jcp('i')|<}=o, because for large d the 
electrons near g = completely dominate the transmis- 
sion. The negative value of Vcp in Fig. [2l^a) is caused by 
the negative polarization {pi < 0) at low barrier height 
Uh for electrons with small q. Vcp remains finite for in- 
finitely high or wide barrier, however, the time to build 
up this voltage, the RC time (tag), goes to infinity due 
to the exponential growth of the resistance. 



When mi is perpendicular to the precession axis of 
m2, i.e. X = 90°, the charge pumping voltage oscil- 
lates around zero because both dot products in Eq. (j7b[) , 
mi • (m2 X 1x12) « 27r/sin6'cos(27r/i) and mi • rh2 = 
27r/sin6'sin(27r/<), give rise to an AC signal. With 
Vcp = a mi • (m2 x 1x12) + b mi • m2, where the two 
components are out of phase by 7r/2, the amplitude is 
given by V^p = 2TTf sinO^/a^ + . An FMR with 6 = 5° 
and / = 20 GHz then gives an AC pumping voltage with 
amplitude as large as V^cp — 0.25/iV. Fig. ^h) shows the 
barrier height dependence of amplitude Vcp quite similar 
to the DC case in Fig.[2)Ja) and for similar reasons. For a 
half metallic junction, the magnitude of the DC pumping 
voltage Vcp can be shown to be bounded by {huj/2e) sin^ 6 
and AC pumping voltage amplitude V^p must be smaller 
than {hiu/2e) shi9, where lo = 27r/. 

When an N layer of thickness a is inserted, it is inter- 
esting to inspect the two different modes, mode 1: FNIF 
and mode 2: FNIF, where F indicates the F layer under 
FMR. Eq. (j7bp applies to mode 1, and applies to mode 
2 with subscript 1 and 2 swapped. The N layer forms 
a quantum well for spin-down electrons that causes the 
oscillation in the charge pumping voltage as a function 
of a as shown in Fig. [3) The period of the quantum oscil- 
lation due to the N insertion layer is about Tr/fcp ~ 3 A. 
However, due to the aliasing effect caused by the discrete 
thickness of the N layer,— the observed period should be 
Tr/jfep — 7r/A|, where A is the thickness of a monolayer. 
In mode 1, the quantum well formed by the N layer can 
modulate T^{a) such that the electrons contributing to 
the transmission the most have Tj~ > (pi > 0) or 
Tj" < (pi < 0), and thus change the sign of the pump- 
ing voltage Vcp. On the other hand, there is no sign 
change in mode 2 because is independent of a. Sim- 
ilar oscillations could also be found for the amplitude of 
the AC pumping voltage. 

Because the AC voltage is proportional to sin 6* and 
DC voltage is proportional to sin^ 9, the AC pumping 
voltage is much larger than the DC counterpart at small 
9. However, in order to observe an AC pumping voltage, 
the time to build up the voltage, the RC time trc = RC, 
has to be shorter than the pumping period, i.e. trc < 
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1//. Approximately trc - (£eo/i/2e2)eV2'"bi/b/?«^d/rf^ 
where e and eo are the dielectric constant and electric 
constant, respectively. A more accurate estimation of 
the RC time for a typical structure is as follows: the 
resistance-area (RA) value of the MTJ in our calcula- 
tion is i?A « 3 n^im^ for d = 0.8 nm (RA « 70 flfim^ 
for d = 1.2 nm, which is consistent with experimental 
values,^) The capacitance of an MgO tunnel barrier with 
d = 0.8 nm is calculated by C/A — eeo/d w 0.1 F/m^ 
(e w 9.7 for MgO). Therefore trc = {RA){C/A) w 0.3 ps 
^/ f ^ 10^ ps. The electromagnetic response is there- 
fore sufficiently fast to follow the AC pumping signal. 

We ignored interface roughness and barrier disor- 
der in the calculation of the pumping voltage. This 
may be justified by the high quality of epitaxial MgO 
tunnel barrier.'^''* Furthermore, the geometric interface 
roughness mainly reduces the nominal thickness of the 
barrier^'* which can be taken care of by an effective thick- 
ness parameter. Impurity states in the barrier open addi- 
tional tunneling channels with [/^ < Uh, which generally 
increases tunneling but also reduces the spin-dependent 
effects when spin-fiip is involved. In general, interface 
roughness and disorder can be important quantitatively, 
but have been shown not to qualitatively affect the fea- 
tures predicted by a ballistic modeli^ In order to be 



quantitatively reliable, the real electronic structure has 
to be taken into account as well. Both band struc- 
ture and disorder effects can be taken into account by 
first-principles electronic structure calculation as demon- 
strated in for metallic structures. 

Recently, a magnetization-induced electrical voltage of 
the order of /iV was measured for an FIN structure by 
Moriyama et alM- The authors explain their findings by 
spin pumping, but note that the signal is larger than 
expected. An FMR generated electric voltage genera- 
tion up to lOO^V was theoretically predicted for such 
FIN structuresi^ Surprisingly, this voltage is much larger 
than hu}/2e ~ /LtV, the maximum "intrinsic "energy scale 
in spin-pumping theory. 

To summarize, a scattering matrix theory is used to 
calculate the charge pumping voltage for a magnetic mul- 
tilayer structure. An experimentally accessible charge 
pumping voltage is found for an FIF MTJ, the pumping 
voltage can be either DC or AC depending on the mag- 
netization configurations. In FNIF structure we find on 
top of the previously reported oscillating TMR^ a charge 
pumping voltage that oscillates and may change sign with 
the N layer thickness. 

This work has been supported by EC Contract IST- 
033749 "DynaMax" . 
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